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ABSTRACT 

We compute the ionizing radiation field at low redshift, arising from Seyferts, 
QSOs, and starburst galaxies. This calculation combines recent Seyfert luminosity 
functions, extrapolated ultraviolet fluxes from our IUE-AGN database, and a new 
intergalactic opacity model based on Hubble Space Telescope and Keck Lya absorber 
surveys. At z = for AGN only, our best estimate for the specific intensity at 1 Ryd 
is Io = l-3lo5 x er g s cm -2 s -1 Hz -1 sr _1 , independent of Hq, Qq, and A. 

The one-sided ionizing photon flux is $i on « 3400± 2 ^ photons cm s , and the 
H I photoionization rate is Thi = 3.2^ 2 '2 x 10 -14 s -1 , for a s = 1.8. We also derive 
Thi for z = — 4. These error ranges reflect uncertainties in the spectral indexes 
for the ionizing EUV (a s = 1.8 ± 0.3) and the optical/UV (auv = 0.86 ± 0.05), the 
IGM opacity model, the range of Seyfert luminosities (0.001L* — 100L*), and the 
completeness of the luminosity functions. Our estimate is a factor of three lower than 
the most stringent upper limits on the ionizing background (3>ion < 10 4 photons cm~ 2 
s _1 ) obtained from Ha observations in external clouds, and it lies within the range 
implied by other indirect measures. Starburst galaxies with a sufficiently large Lyman 
continuum escape fraction, (/ esc ) > 0.05, may provide a comparable background to 
AGN, Iq(z = 0) = l-l^o 7 x 10~ 23 ergs cm -2 s _1 Hz^ 1 sr" 1 . An additional component 
of the ionizing background of this magnitude would violate neither upper limits from 
Ha observations nor the acceptable range from other measurements. 



Subject headings: intergalactic medium — diffuse radiation — galaxies: Seyfert - 
galaxies: starburst 
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1. INTRODUCTION 

The ionizing background that permeates intergalactic space is of fundamental interest for 
interpreting QSO absorption lines and interstellar high-latitude clouds. Produced primarily by- 
quasars, Seyfert galaxies, and other active galactic nuclei (AGN), these Lyman-continuum (LyC) 
photons photoionize the intergalactic medium (IGM), set the neutral hydrogen fraction in the Lya 
forest absorbers, and help to determine the ion ratios in metal-line absorbers in QSO spectra. 
Ionizing radiation may control the rate of evolution of the Lya absorption lines at z < 2 (Theuns, 
Leonard, & Efstathiou 1998; Dave et al. 1999), and it may affect the formation rate of dwarf 
galaxies (Efstathiou 1992; Quinn, Katz, & Efstathiou 1996). The hydrogen photoionization rate, 
Thi(^), is an important component of N-body hydrodynamic modeling of the IGM. Because of the 
large photoionization corrections to the observed H I absorption, the inferred baryon density of 
the IGM and metal abundance ratios also depend on the intensity and spectrum of this radiation. 
Within the Milky Way halo, the ionizing background can affect the ionization state of high-velocity 
clouds located far from sources of stellar radiation (Bland-Hawthorn & Maloney 1999). 

The ionizing background intensity at the hydrogen ionization edge {hvQ = 13.6 eV) is denoted 
Iq, in ergs cm~ 2 s _1 Hz -1 sr _1 , hereafter denoted "UV units" or understood in context. In 
an optically thin environment, the background spectrum reflects that of the sources, QSOs and 
Seyfert galaxies, which appear to have steep EUV spectra of the form F v oc (i//vo)~ as , with 
a s = 1.77 ± 0.15 from 350-1050 A (Zheng et al. 1997), or starburst galaxies with a s ss 1.9 - 2.2 
(Sutherland & Shull 1999). At high redshift, the IGM is optically thick, owing to the numerous 
Lya absorbers that ionizing photons must traverse. Thus, the background spectrum is strongly 
modified by absorption and re-emission (Haardt & Madau 1996; Fardal, Giroux, & Shull 1998, 
henceforth FGS). 

Estimates of Iq at high redshift are usually obtained from the "proximity effect" (Bajtlik, 
Duncan, & Ostriker 1988), the observed paucity of Lya absorbers near the QSO emission redshift. 
Recent measurements give values of Iq ss 10~ 21 UV units at z ss 3: Io = I.OIqJ x 10~ 21 (Cooke, 
Espey, & Carswell 1997), I = (0.5 ± 0.1) x 10~ 21 (Giallongo et al. 1996), and I = 0.75 x 10~ 21 
(Scott et al. 1999). At low redshift, the lower comoving density of QSOs and their diminished 
characteristic luminosities suggest that the metagalactic background is reduced by about a 
factor of 10 2 to io ss 10~ 23 . Using an optical QSO luminosity function (cf. Boyle 1993) and an 
empirical model of IGM opacity (Miralda-Escude h Ostriker 1990), Madau (1992) estimated that 
Io = 6 x 10~ 24 at z = 0. However, this is an uncertain estimate, which now appears low compared 
with several local determinations. Theoretical extrapolations of Iq to low z are uncertain because 
they depend sensitively on the assumed AGN luminosity function and on the IGM opacity model 
(Giallongo, Fontana, & Madau 1997; FGS). As we will show, the low-z IGM opacity appears to 
be dominated by Lya absorbers in the range 14 < log Nhi < 18, for which we are just starting 
to obtain statistically reliable information from the Hubble Space Telescope (HST). A future key 
project on Lya absorbers with HST and a Ly/3 survey with the Far Ultraviolet Spectroscopic 
Explorer (FUSE) should be even more enlightening. 
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Observational estimates of or upper limits on Iq at low redshift have been made by a variety 
of techniques, as described in Table 1. These methods include studies of the proximity effect at 
(z) ~ 0.5 (Kulkarni & Fall 1993), edges of H I (21 cm) emission in disk galaxies (Maloney 1993; 
Dove & Shull 1994a), and limits on Ha emission from high-latitude Galactic clouds (Vogel et al. 
1995; Tufte, Reynolds, & Haffner 1998) and extragalactic H I clouds (Stocke et al. 1991; Donahue, 
Aldering, & Stocke 1995). Since all these techniques are based on the integrated flux of LyC 
radiation, it is convenient to define <I>ion (photons cm~ 2 s _1 ), the normally incident photon flux 
through one side of a plane. For an isotropic, power-law intensity, I v = Iq{u / Vo)~ as , we can relate 
the integral quantity <&i on to the specific intensity Io : 

<D ion = 2*j\d» £° ^ dv = (J^j = (2630 cm- 2 s" 1 )/^ (^) , (1) 

where n = cos is the angle relative to the cloud normal and i_23 is the value of Iq expressed in 
units of 10~ 23 UV units. Most of the upper limits on Iq translate into values of <& lon in the range 
10 4 - 10 5 photons cm s . For an assumed EUV spectral index a s ~ 1.8, and the approximate 
form, o v ~ ao(v / vo)~ 3 , for the H I photoionization cross section, the hydrogen photoionization 
rate due to this metagalactic intensity is, 

r ™ " = < 2 - 49 x 10 ~ 14 s_1 ) / - (^r-) • < 2) 

h(3 + a s ) \3 + a s J 

where do = 6.3 x 10~ 18 cm 2 is the hydrogen photoionization cross section at hv®. 

These Ha measurements and limits are improving with better Fabry-Perot techniques 
(Bland-Hawthorn et al. 1994; Tufte et al. 1998). In addition, we now have more reliable HST 
measurements of the opacity from the low-redshift Lya clouds (Weymann et al. 1998; Shull 1997; 
Penton, Stocke, & Shull 1999). Therefore, better computations of the metagalactic radiation field 
are timely. In this paper, we compute the contribution of Seyfert galaxies, QSOs, and starburst 
galaxies to the low-redshift ionizing background using three ingredients: (1) a Seyfert/QSO 
luminosity function; (2) AGN fluxes at A < 912 A from extrapolated IUE spectra; and (3) an 
improved IGM opacity model, based on recent HST surveys of Lya clouds at low redshift. In 
§ 2 we describe these ingredients. In § 3 we give the results for Iq and § lon at z ~ 0, together 
with error estimates. In § 4 we summarize our results and discuss future work that could improve 
estimates of Iq. 



2. METHODOLOGY 

The solution to the cosmological radiative transfer equation (Peebles 1993) for sources with 
proper specific volume emissivity e(u, z) (in ergs cm -3 s _1 Hz -1 ) yields the familiar expression 
(Bechtold et al. 1987) for the mean specific intensity at observed frequency v b s as seen by an 
observer at redshift z Q b s : 

Iv{v bs,z ohs ) = — I e(u,z) exp(-T eff ) dz . (3) 

47r Jz obs az [L + z) 
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rz roo Q 2 J\f 

Te&(vobs,z bs,z) = dz'l — — [I - exp(-r)] dN HI , (4) 



Here, v = f D fes(l + + z G &s) is the frequency of the emitted photon (redshift z) observed at 
frequency (redshift z b s ), dl/dz = (c/Hq)(1 + z)~ 2 (l + VLqz)~ 1 I 2 is the line element for a 
Friedmann cosmology, and r c g is the effective photoelectric optical depth due to an ensemble of 
Lya absorption systems. For Poisson-distributed clouds (Paresce, McKee, &; Bowyer 1980), 

dN HI dz' 

where d 2 N /dNujdz' is the bivariate distribution of Lya absorbers in column density and redshift, 
and r = Nhi cr{u) is the photoelectric (LyC) optical depth at frequency u due to H, He I, and 
He II through an individual absorber with column density Nhi- For purposes of assessing the 
local attenuation length, it is useful (Fardal & Shull 1993) to use the differential form of eq. (4), 
marking the rate of change of optical depth with redshift, 

dr eS (°° 8 2 M 

-dz- = L dN^dz [l - QM - T)]dNHI - (5) 

The attenuation length, in redshift units, is given by the reciprocal of dr e g/dz. At low z, since 
dr c fi/dz <J 1 at the hydrogen threshold, its frequency dependence is significant, and the attenuation 
length can extend to z ~ 2. In the past few years, more sophisticated solutions to the cosmological 
transfer have been developed (Haardt & Madau 1996; FGS) taking into account cloud emission 
and self-shielding. Figure 1 illustrates our group's recent calculation of the ionizing background 
spectrum, computed in full cosmological radiative transfer with a new IGM opacity model based 
on high-resolution Keck spectra of the Lya forest and local continua. These models include cloud 
self-shielding and emission. We have connected this high-redshift opacity model with our new 
model from HST studies (discussed below) for the low-redshift opacity at a transition redshift 
z = 1.9. By redshift z = 0, the intensity has declined to Iq ps 1.3 x 10~ 23 , corresponding to 
$i on ss 3000 photons cm~ 2 s for sources with 

In the work that follows, we compute I v using both our detailed cosmological radiative 
transfer code and an approximate solution to equation (3). In this approximation, we neglect 
the effects of emission from attenuating absorbers and approximate the opacity with a simple 
power-law fit that neglects the effects of He absorption. We will discuss the accuracy of this 
approximation in more detail in § 3. Because the opacity of the IGM is much smaller at low 
redshift, this more rapid calculation is adequate for estimating the present-day level of radiation 
just above uq. The primary ingredients for the computation of Iq at low redshift are the source 
emissivity, e(i/,z), and the opacity model for T e ff (v b s , z bs, z )- I n the following sub-sections, we 
describe how we determined these quantities. 



2.1. AGN Luminosity Function and Spectra 



The distribution of AGN luminosities is typically described by a rest-frame, B-band luminosity 
function. In order to estimate the total emissivity of AGN at the Lyman limit, we must know 
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both the luminosity function and the average spectrum of the AGN. In addition, we must know 
the assumptions about the spectrum that were made to construct the luminosity function. 

To estimate the intrinsic AGN quasar spectrum, we begin with the Seyfert optical sample 
of Cheng et al. (1985), based on Seyfert 1 and 1.5 galaxies covered by the first nine Markarian 
lists. Their sample was corrected for incompleteness, and the contribution from the host galaxy 
was subtracted out. The separation of nuclear and host galaxy luminosity becomes increasingly 
challenging at the faint end of the luminosity function. Ideally, careful, small-aperture photometry 
would be used for these estimates. Cheng et al. (1985) relied instead on two independent methods 
to separate the contribution from the nucleus. In the first method, they assumed a template host 
galaxy and corrected this for orientation and internal extinction. In the second method, they 
assumed that all nuclei had the same intrinsic colors and determined the nuclear contribution 
via the color-given method of Sandage (1973). They found these two methods to give consistent 
results. In addition, they compared the color-given method with nuclear magnitudes derived by 
careful surface photometry for 11 Seyfert 1 galaxies (Yee 1983). They assigned a total uncertainty 
of 0.5 mag in the nuclear Mb to the sample. On re-examination of the sample, we found that the 
errors most likely decrease with the luminosity of the galaxy. We assume that the errors on the 
specific luminosity, Lb (ergs s -1 Hz -1 ), decrease linearly from 0.24 dex at log Lb = 28 to 0.16 dex 
at log Lb = 30. 

From this sample of Seyferts, we chose 27 objects observed repeatedly by the International 
Ultraviolet Explorer (IUE) satellite. Together with many other AGN, these Seyfert galaxies are 
part of the Colorado IUE- AGN database (Penton, Shull, & Edelson 1996), which gives both mean 
and median spectra. Since these AGN are subject to flux variability, the distribution in flux is a 
skewed distribution with a tail that includes short flares studied by various IUE campaigns. To 
provide a conservative estimate of the ionizing fluxes, we have therefore used median spectra to 
derive correlations; however the differences in the correlations are only a few percent. The line- free 
regions of the median IUE spectra were fitted to power-law continua and extrapolated to 912 A 
(rest-frame), from which we derive the specific luminosity, L912 (ergs s -1 Hz -1 ). We also convert 
from M B to L B , at v = 6.81 x 10 14 Hz (4400 A) by the formula derived from Weedman (1986, eqs. 
3.15 and 3.16), log Lb = 0.4(51.79 — Mb)- Figure 2 shows the correlation between L912 and Lb- 
The error bars are only shown for Lb, since the errors in L912 are much smaller. We find that 
L 9 i 2 = (2.60 ± 0.22) x 10 28 (L B /10 29 )( L114±0 - 081 ) erg s -1 Hz -1 . This translates to a UV-optical 
spectral slope that depends on Lb as 

a uv = (0.86 ± 0.05) - (0.16 ± 0.12) log(L B /10 29 ) . (6) 

The evidence for Lb dependence is marginal, however, and our basic model will simply assume a 
constant slope auv = 0-86 between the B band and the Lyman limit. This agrees quite well with 
the average QSO spectrum derived by Zheng et al. (1997). 

The second ingredient in the computation of the AGN emissivity is the luminosity function. 
The function that matches observations over the broadest redshift and luminosity range is the 



-6- 



comoving analytic form given by Pei (1995, eqs. 6-8), 



nL ' Z} = {L/L z y l + {L/L z )h ' (7j 
where the characteristic "break luminosity" is given by 

L z = L^l + zJ-^-^exp [-(z-z*) 2 /2oi 

= L (l + z)- (1 - auv) exp[-z(z-2^)/2cj2] . (8) 

Here, «uv is the UV-optical spectral index and Lq = L* exp(— zl/2a1) is the "break luminosity", 
L z , at the present epoch. Note that we define our spectral index with the opposite sign to 
that assumed by Pei, i.e., we adopt L u oc v~~ auv where most data suggest that 0.5 < aijv < 1- 
We present results for luminosity functions based on the two sets of assumptions about the 
cosmology and optical spectral index cktjv derived by Pei. The "open model" has h = 0.5, 
n = 0.2, and a vy = 1.0 and yields 0i = 1.83, h = 3.70, z* = 2.77, log(L*/L©) = 13.42, 
(j* = 0.91, and log(<£*/Mpc~ 3 ) = —6.63. The "closed universe" model has h = 0.5, 1\) = 1, 
and a uv = 0.5 and yields = 1.64, (3 h = 3.52, = 2.75, log(L*/L ) = 13.03, a* = 0.93, and 
log(<E>*/Mpc~ 3 ) = —6.05. Figure 3 shows these two models at z = 0. After correcting for the 
different spectral indices, there is no physical reason why these two models should differ in the 
ionizing intensity they imply or in their value as z — > 0. However, the intensity estimates from 
these models differ by up to 40% (see FGS). This points to substantial uncertainties in the Pei fit 
that are not reflected in the small formal errors. 



max • 



We integrate the ionizing luminosity density over the luminosity function from L rmn to L 
In our standard model, we assume that these limits are O.OIL^ and WL Z , respectively, and we 
explore the sensitivity of the results to L m \ n . If we define x = L/L z , we can write the comoving 
specific volume emissivity as 



where we assume an EUV power-law spectral index a s for the AGN and integrate from 
Xmin = L m \ n /L z to x max = L max / L z . With x m j n = 0.01 and x max = 10, our results are insensitive 
to x max and moderately insensitive to x m \ n . An increase (decrease) of x imn by a factor of 3 
leads to a decrease (increase) of the calculated emissivity by 15%. A similar change of x max 
changes the emissivity by only 2%. From several trials, we find an empirical scaling relation, 
Iq^ gn oc (L m i n /0.01L*)~°' 17 (a s /1.8)~ ' 97 . To better assess the uncertainty from the extrapolation 
to bright and faint sources, we use Figure 2 of Pei (1995), which shows the range of B luminosities 
that contribute to the luminosity function fit. We define the "completeness" as the integral of the 
emissivity over this range, compared with the integral for the standard range 0.01 < x < 10. This 
completeness rises from 20% at z = to 80% at z ~ 2, but falls to 20% at z ~ 4. The low-z results 
can be made more robust by considering additional surveys of Seyferts and QSOs. 
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Figure 3 shows the B-magnitude luminosity function from the Cheng et al. (1985) sample, 
as well as the results from the survey of Kohler et al. (1997) for Seyferts and QSOs with z < 3. 
These results are also compatible with the results of the local optical luminosity function based on 
an X-ray selected sample of AGN (Delia Ceca et al. 1996). There are some discrepancies, however, 
between these results and the Pei luminosity function. There are more QSOs at the bright end of 
the luminosity function at z < 1 than in the Pei form, an effect that has been noted by several 
authors (Goldschmidt et al. 1992; Goldschmidt k Miller 1998; La Franca & Cristiani 1997). This 
appears to have its origin in systematic errors in the Schmidt & Green (1983) QSO survey. In 
addition, there appears to be a slight deficit at the knee of the Pei function in both the Kohler et 
al. and Cheng et al. samples. In fact, the Kohler results are fitted adequately with a single power 
law, as shown in Fig. 3. 

It turns out that these two effects cancel to within 5% when we compute the total B-band 
emissivity over the range 0.01 < Lb/L z (0) < 10. However, the fact that the luminosity function 
appears to be changing shape weakens the assumptions leading to the Pei function, and the 
convergence at the faint end becomes even slower. 

Hence, our best estimate of the emissivity is unchanged from the Pei (1995) model, but we 
suspect that there are still substantial uncertainties in the AGN luminosity function. In Fig. 4 we 
plot the emissivity from the Pei "open" model. To find a reasonable range of emissivity models, 
we multiply (or divide) this by the square root of the "completeness" of the Pei model as derived 
above, and we also change the spectral slope cktjv and its dependence on Lb by one standard 
deviation to minimize (or maximize) the ionizing emissivity. We consider the band thus defined to 
be a conservative la range for the emissivity. Note that the Pei "closed" model, adjusted for the 
different spectral index and cosmology, lies well within this range. 

The average UV QSO spectrum derived by Zheng et al. (1997) implies an EUV spectral 
index a s = 1.77 ± 0.15 for radio-quiet QSOs and a s = 2.16 ± 0.15 for radio-loud quasars. Since 
radio-quiet QSOs are much more numerous, we select a s = 1.8 ± 0.15 as a representative measure 
of the EUV spectral index. 

2.2. Stellar Contributions to the Ionizing Emissivity 

The contribution of stars within galaxies to the ionizing background remains almost 
completely unknown at all epochs. It has been realized (cf. Madau & Shull 1996) that the amount 
of ionizing photons associated with the production via supernovae of the observed amount of 
metals in the universe might easily exceed the ionizing photons produced by AGNs. The problem 
has been to estimate the average fraction, (f eS c), of ionizing photons produced by O and B stars 
that escape the galaxies into the IGM. Thus far, observational limits on (f esc ) exist from only a 
few nearby galaxies (Leitherer et al. 1995; Hurwitz et al. 1997). On theoretical grounds, Dove & 
Shull (1994b) concluded that an escape fraction of order 10% might be possible, while more recent 
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models of the escape of ionizing photons through supershell chimneys (Dove, Shull, & Ferrara 
1999) suggest fractions of 3-6%, which are compatible with the observational limits. Recently, 
Bland-Hawthorn & Maloney (1999) used measurements of Ha from the Magellanic Stream to infer 
an escape fraction ~ 6% from the Milky Way. 

Deharveng et al. (1997) argued that even the present indirect upper limits on Iq must limit 
the present-day escape fraction to below 1%. They depart radically, however, from the work 
of Dove & Shull (1994b) and Dove et al. (1999) in their treatment of the nature of the escape 
fraction. Deharveng et al. assume that stellar ionizing photons are prevented from escaping their 
host galaxies by the opacity of neutral hydrogen, effectively in the limit that the hydrogen forms 
an unbroken sheet. Consequently, the opacity varies approximately as (f/i/o) -3 , so that the escape 
of higher-energy photons is dramatically enhanced. For example, the optical depth can decrease 
from r ~ 50 at the Lyman edge to r < 1 at 4 rydbergs. For photons emitted at high redshift that 
survive to contribute to Iq at z = 0, the sharp increase in the escape fraction outweighs the effects 
of redshifting. In the Deharveng picture, the contribution of high-z galaxies to the present-day 
mean intensity at v = v$ is strongly enhanced. Thus, if (/ esc ) ~ 0.001, starbursts may match the 
contribution of quasars to the present-day ionizing mean intensity. 

However, we believe the Deharveng model for photon escape is physically unrealistic. Our 
alternative view is that the internal galactic opacity to all stellar photons is large, and photons 
may only escape from isolated regions of high star formation whose H II regions or attendant 
supershells are able to break through this high opacity layer (Dove et al. 1999). Within this view, 
a constant escape fraction with frequency is more appropriate. 

Our estimate of the stellar ionizing photons is made in the following way. Gallego et al. (1995) 
performed an Ha survey of galaxies and fitted their derived luminosity function to a Schechter 
function. By integrating this, they were able to estimate a total Ha luminosity per unit volume 
at low redshift, Lhu = 10 3ai±0 ' 2 erg s _1 Mpc~ 3 . In the usual fashion, we relate this to the 
number of ionizing photons by dividing the number of Ha photons by the fraction of the total H 
recombinations that produce an Ha photon. We multiply their estimate of the ionizing photon 
production by (/ esc ) to give the total ionizing photons in the IGM. The representative models 
assume (/ esc ) = 0.05 and are shown in Figure 5. Preliminary results are now available from 
the KPNO International Spectroscopic Survey (KISS) (Gronwall et al. 1998), which probes to 
fainter magnitudes than the Gallego et al. survey. Gronwall (1998) quotes a value Lu a = 10 39 ' 03 
erg s" 1 Mpc -3 , but notes that these results are preliminary and represent a lower limit to the 
true Ha density, since even their deeper survey is incomplete for galaxies with Ha emission-line 
equivalent widths less than 25 A. 

The spectrum of ionizing photons from clusters of hot stars differs from that of AGN in that 
stars emit relatively few photons more energetic than 4 Ryd. Sutherland & Shull (1999) have 
shown that, between 1 and 4 Ryd, the spectrum of a starburst may be approximated as a power 
law with spectral index a s = 1.9 — 2.2. 
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We have tied the redshift evolution in ionizing photon production rate to the star formation 
evolution observations of Connolly et al. (1997) based upon the Hubble Deep Field (HDF). The 
effects of dust extinction remain a major uncertainty in the determination of the star formation 
rate at high redshift (cf. Pettini et al. 1998). Many high-z star- forming galaxies appear to be 
dust-obscured, based on recent sub-millimeter studies of the HDF (Hughes et al. 1998; Barger 
et al. 1998). In addition, a survey of Lyman-break galaxies at z 3 — 4 (Steidel et al. 1999), 
covering a much larger angular extent than the HDF, finds no significant difference in the star 
formation rate at z = 3 and z = 4. With corrections for dust extinction, the star formation could 
remain constant from z = 1.5 out to z > 4. As a result, we have also considered a case in which 
the star formation remains constant after reaching its peak at z ~ 2. However, despite the many 
uncertainties of high-z star formation rates, the effects on the present-day level of ionization are 
minimal, of order a few percent. 



2.3. Absorption Model for the IGM 

At z < 2, the redshift densities of the Lya clouds and Lyman limit systems decline steeply 
with cosmic time. Morris et al. (1991) and Bahcall et al. (1991) used HST observations to show 
that this decline could not be extrapolated to the present, as far too many Lya absorbers were 
observed toward 3C 273. The large data set gathered by the Hubble Key Project with the Faint 
Object Spectrograph (FOS) shows a sharp break in the evolving redshift density at z ~ 1.5 — 2 
(Weymann et al. 1998). Ikeuchi & Turner (1991) showed that the cessation in this steep decline 
was a natural consequence of the falloff in the ionizing emissivity from z = 2 down to 0. This 
conclusion has been borne out by detailed cosmological simulations (Dave et al. 1999), which 
indicate that the effect is insensitive to the specific cosmological model. 

In our calculations, we will base our absorption model on observations. Our analysis follows 
the traditional "line-counting" method, where spectral lines are identified by Voigt profile-fitting 
and the opacity is calculated by assuming a Poisson distribution of these lines. At high redshift, 
FGS used this method to estimate the opacity based on high-resolution observations of QSO Lya 
absorption lines from Keck and other large-aperture telescopes in the redshift range 2 < z < 4. It 
is not appropriate to extrapolate this model to lower redshifts, owing to the rapid evolution rate 
of the Lya forest. 

Our method of determining dr e s/dz considers Lya lines in three ranges of column density: 
from 12.5 < log Nhi < 14.0 (HST/GHRS survey), from 14.0 < logN HI < 16 (HST/FOS survey), 
and for logNjji > 17 (HST/FOS Lyman-limit survey). The HST/FOS survey forms the core 
of our standard opacity model. We combine these results with HST/GHRS measurements of 
weak lines and with the HST/FOS survey of Lyman-limit systems, extrapolated downward from 
the Lyman limit by two different methods. At redshifts z < 1.5, the most extensive study of 
strong Lya absorbers (10 14 ~ 16 cm~ 2 ) is the QSO Absorption Line Key Project with HST/FOS 
(Januzzi et al. 1998; Weymann et al. 1998). Weaker Lya lines, which contribute a small amount 
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to the opacity and serve as a constraint on the column density distribution, were studied by Shull 
(1997), Shull, Penton, & Stocke (1999), and Penton et al. (1999) using HST/GHRS spectra. The 
distribution of Lyman limit systems with Nhi > 10 17 cm~ 2 are discussed by Stengler-Larrea et 
al. (1995) and Storrie-Lombardi et al. (1994). Each of these surveys suffers from incompleteness 
or saturation effects in various regimes. Therefore, extrapolations outside the range of Nhi and 
comparisons in regimes of overlap are helpful. Extensive future UV surveys with HST and FUSE 
will also reduce some of the uncertainties (see discussion in § 4) . 

HST/GHRS studies of Lya absorbers in the range 12.5 < logN HI < 14.0 (Penton et al. 
(1999) find a column density distribution, dM /dNni oc A^}' 74±0 ' 26 . The cumulative opacity of 
these weak lines (up to 10 14 cm~ 2 ) is relatively small: dr c fi/dz = 0.025 ± 0.005 for the low-redshift 
range (0.003 < z < 0.07). However, a small number of higher column density systems produce a 
steady rise in the cumulative opacity for log Nhi > 14. Extending the HST/GHRS distribution 
up to log Nhi = 15 gives dr e s/dz ss 0.09 ± 0.02. Above this column density, saturation effects and 
small-number statistics make the number counts more imprecise. In the range 15 < log Nhi < 16, 
Penton et al. (1999) estimate an additional contribution, dT c g/dz 0.1 — 0.3. 

The HST/FOS Key Project spectra have insufficient resolution to determine line widths or to 
resolve velocity components. As a result, the conversion from equivalent width, W\, to column 
density, Nhi, is difficult for saturated lines. In the absence of other lines (e.g., Ly/3), one can only 
estimate the conversion from W\ to Nhi by assuming a curve of growth and doppler parameter. 
This difficulty was noted by Hurwitz et al. (1998) in their comparison of unexpectedly strong 
ORFEUS Ly/3 absorption compared to predictions from HST Lya lines toward 3C 273. For 
unsaturated lines, W\ = (54.4 mA)A"i3, where Nhi = (10 13 cm _2 )A / "i3 and where the line-center 
optical depth is r = (0.303)Aq 3 for a doppler parameter b = (25 km s _1 )6 2 5- The HST/FOS 
Key Project lines with W\ > 240 mA are highly saturated in the range (14 < log Nhi < 17) that 
dominates the Lya forest's contribution to the continuum opacity. 

As a first attempt to incorporate the Key Project information, we focus on the statistical 
frequency of Lya absorbers, dAf/dz = 30.7 ± 4.2, for lines with W\ > 240 mA. For b = 25 km s _1 , 
this corresponds roughly to log Nhi > 14. We compute opacities based upon sample 5 of Weymann 
et al. (1998), which includes 465 absorption lines (W\ > 240 mA) that could not be matched 
with corresponding metal lines. This sample was intended to remove high column density lines 
that may evolve more rapidly with redshift, consistent with the results of Stengler-Larrea et al. 
(1995) for Lyman limit systems. This segregation has little effect on the opacity. The Key Project 
became incomplete at column densities well below the Lyman limit. Therefore, to determine an 
IGM opacity, we needed to extrapolate to log Nhi = 17 using assumptions about the distribution. 

The bivariate distribution of Lya absorbers per unit redshift and unit column density can be 
expressed as d 2 N ' jdz dNni = A(Nhi/W 17 cm~ 2 )~^(l + z) 7 . Improvements on this form have been 
suggested, notably the addition of one or more breaks in the power-law distribution (Petitjean 
et al. 1993; FGS). Here, we parameterize our uncertainty by assuming just one power-law index, 
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but we vary the upper limit on column density to which we integrate. Figure 6 shows a set of 
curves, corresponding to (3 = 1.5 ± 0.2, of the differential effective opacity, dr e s/dz, evaluated 
at z = and at the hydrogen threshold. Assuming (3 = 1.5 and no break in the distribution, 
we find dr e s/dz « 0.2 for the FOS range 10 14 ~ 16 cm~ 2 and dr e s/dz « 0.5 for the expanded 
range 10 14 ~ 17 cm~ 2 . Because 7 = 0.15 ± 0.23 for this FOS sample, the opacity does not change 
greatly with redshift. We choose a standard value dr e ff/dz ss 0.5, corresponding to (3 = 1.5 and 
A max = 10 17 cm" 2 . 



The Lyman-limit data (Stengler-Larrea et al. 1995) can be fitted to the form 

<,TvU 0.263(1 + z) L50 , (10) 



dz 

assuming that d 2 N ' /dzdNni oc A^}' 5 ^ + z) 1 - 50 . In absorption model 1, we took the lower limit 
on Nhi for partial LL systems to be N\ = 10 17 cm~ 2 . If this limit is extended down to 10 16 cm -2 
or 10 15 - 5 cm -2 , the coefficient 0.263 in eq. (10) increases to 0.382 and 0.411, respectively; the 
latter choice becomes our Model 2. 

We summarize our three opacity models in Table 2. We use equation (5) to calculate the 
opacity, but include an approximation of the frequency dependence of the opacity over the range 
1-3 Ryd, which dominates the H I photoionization rate. This approximation takes the form 

dr cS 

dz ' V^o, 

In Figure 7, we compare dT e s/dz at v = v$ for the three models described above, as well as 
low-redshift extrapolations of the opacities of FGS and Haardt & Madau (1996). It can be seen 
that the poorly determined column-density distribution of the Lya forest leaves a large uncertainty 
in the total opacity, even though the evolution of the number density is tightly constrained by the 
HST observations. The partial LL systems (16 < log Nhi < 17.5) probably dominate the IGM 
opacity at low redshift, but they are so rare that statistical fluctuations from sightline to sightline 
are quite large. It will require many high signal-to-noise spectra along low redshift lines of sight 
to reduce the uncertainty in their contribution to the opacity. 



cA-\ (1 + *)-*. (ii) 



3. RESULTS: THE IONIZING RADIATION FIELD 

3.1. The Contribution from AGN 

Our best-estimate model for the present day intensity Iq makes the following four assumptions: 
(1) The AGN distribution is described by the Pei (1995) QSO luminosity function with h = 0.5 
and f2o = 0.2, modified by the assumption that the optical-UV spectral index a\jy = 0.86, but 
with no correction for intervening dust. (2) The lower (upper) cutoffs to the luminosity function 
are £ m i n (max) = 0.01(10)L*. (3) The opacity below z « 1.9 is our "standard model," itself based 
upon HST observations of Lya forest and Lyman Limit systems. Above z « 1.9, it is Model A2 
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from FGS. (4) The ionizing spectrum (y > uq) has spectral index a s = 1.8. Using the full radiative 
transfer calculation outlined in FGS, we find 

I = l.3+g : | x 10~ 23 erg cm" 2 s' 1 Hz" 1 sr" 1 . (12) 

As we discuss below, most of the uncertainties in our estimate are essentially multiplicative. The 
uncertainties quoted above represent the addition in quadrature of the uncertainties in the log of 
the factors discussed below. Direct addition in quadrature of the relative uncertainties gives a 
similar range of uncertainty. Strictly speaking, although none of the uncertainties is independent, 
we treat them as if they were in computing the total uncertainty. 

Figure 4 summarizes the uncertainties in the emissivity due to AGN that pertain to 
assumptions (1) and (2). The corresponding 1 a uncertainty in the calculated specific intensity for 
the £Iq = 0.2 case is ±0.19 dex. As FGS noted, the intensity should be independent of the choice 
of Hq or Qq if the emissivity is actually determined observationally. The AGN luminosity function 
is propagated from low to high redshift, assuming pure luminosity evolution, a parameterization 
that depends only on z, and not explicitly on appropriate cosmological distances and luminosity 
evolution parameters. The difference in the calculated mean intensity for a closed universe 
(Qq = 1) is not an independent source of uncertainty. However, we estimate an uncertainty of 
±0.01 dex, which is negligible. 

The uncertainties inherent in our choice of opacity model appear both in the degree of 
attenuation of ionizing photons due to the absorbers, and in the contribution of diffuse ionizing 
radiation from the absorbers (mainly He II ionizing photons reprocessed to He II Lya and 
two-photon radiation) . Using Models 1 and 2 instead of our standard model, we calculate slightly 
higher levels of Iq, by factors of 0.0068 and 0.0329 dex. A sample standard deviation of the 
models is ±0.017 dex and is probably an adequate assessment of the la uncertainty in Iq due to 
the opacity. A more complete picture of the full (at least several a) uncertainty may be obtained 
from the following extreme cases. The diffuse ionizing radiation from absorbers contributes 20% 
of the ionizing intensity at z = (see Fig. 8a), so eliminating this contribution reduces Iq by 0.1 
dex. If there is no opacity at all, Iq increases by 0.26 dex. If the column density distribution for 
the absorbers in our standard model is assumed have an unbroken power law with (3 = 1.3 rather 
than 1.5, the reduction in Iq is 0.21 dex. The shape of the spectrum of AGN shortward of 912 A 
is a relatively small source of uncertainty. As a s is varied between 1.5 and 2.1, representative of 
the 2 a uncertainty in a s , the relative specific intensity varies by ±0.1 about the standard value 
a s = 1.8. We therefore adopt a conservative uncertainty of ±0.1 dex in Jo- 
There remains an additional systematic uncertainty in the contribution of AGN to Iq. As 
discussed by FGS, the Pei luminosity function used here produces insufficient ionizing photons 
to account for the level of the ionizing background at z > 3.5 implied by the proximity effect. 
Further, there are not enough ionizing photons to reionize the universe by z ~ 5 (Madau 1998), 
the epoch of the highest redshift quasars. While the latter point remains a problem for scenarios 
in which only AGN ionize the IGM, the former difficulty is ameliorated by the suggestion that, 
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especially at higher redshifts, the number of AGN are being unelercounted due to obscuration by 
dust-laden absorbers (Heisler & Ostriker 1988; Fall & Pei 1993; Pei 1995). Using the dust-corrected 
luminosity function of Pei (1995), we find that 1$ at z = is increased by 0.08 dex. 

As shown in Fig. 1, the level of the mean intensity at z ph 2 calculated by FGS is Iq ph 7x 10~ 22 . 
For this model, ionizing radiation due to sources at z > 2 produces 20% of the mean intensity at 
z = (see Fig. 8a). As Fig. 1 shows, redshifted He II Lya diffuse emission is still substantial 
at z = 2. Because He II A304 emitted at z > 2 is redshifted below the H I threshold by z = 0, 
this emission contributes less than 10% at z = 0. To give a specific example, suppose that the 
number of AGN at high redshift was severely underestimated, so that the metagalactic background 
due to AGN at z ~ 2 — 3 was quintupled, while retaining the same spectral shape. Then, the 
value of Iq at z = would be increased by only 0.15 dex. Even for AGN, the strong attenuation 
of high-energy photons by He II absorption greatly limits any contribution to the present-day 
ionizing background from sources at z > 3. A significantly larger population of AGN at z = 2 — 3 
will not augment Iq by more than 40% at z = 0. This systematic effect has not been included in 
the uncertainty in eq. (12). 

3.2. The Possible Contribution from Hot Stars 

An estimate of the ionizing radiation contributed by hot stars is complicated by its dependence 
on factors such as (/ eS c}j f° r which a good estimate of its magnitude and uncertainty does not 
exist. As a result, we include these factors explicitly in our results. Our model for the possible 
contribution of stars to the present day specific intensity Iq makes the following assumptions: (1) 
The production of ionizing photons by stars at the present time may be calibrated by the density 
of Ha photons in the extragalactic background. (2) The star formation rate is proportional to that 
from the observations of Connolly et al. (1997) and assumes h = 0.5 and fio = 0.2. (3) The IGM 
opacity is our "standard model." (4) The average spectrum of the OB associations that provide 
the ionizing photons has spectral index a s = 1.9 (1-4 Ryd) with no photons above 4 Ryd. Using 
a radiative transfer calculation that neglects the diffuse radiation contributed by intervening 
absorbers (see discussion below), we find 

I = 1.1±- x 10-- (^£>) ( J*«l) erg cm- s" 1 Hz"* sr^ , (13) 

where we have scaled the Ha luminosity density, Ln a , to the Gallego et al. (1995) standard value, 
10 39 1 erg s" 1 Mpc~ 3 , and adopted a probable LyC escape fraction, (/ esc ) = 0.05. 

The uncertainties for which we have not directly parameterized our ignorance are mainly 
mutiplicative, as in the case for AGN. The uncertainty in the numerical coefficient quoted 
above again represents the addition in quadrature of the uncertainties in the log of the factors 
discussed below. We make the same conservative assumption that the individual uncertainties are 
independent. The uncertainty in the Ha luminosity density suggested by Gallego et al. (1995) 
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is ±0.2 dex, while the uncertainty in the preliminary KISS result, is ±0.1 dex. We adopt an 
uncertainty of ±0.2 dex for the Ha calibration, but it is possible that this uncertainty will be 
reduced substantially when full details of KISS are released. As long as (/ esc ) is small, so that 
recombinations within the galaxies provide a fair accounting of the number of ionizing photons 
produced, this is the uncertainty we associate with assumption (1). The conversion of Ha photons 
to ionizing photons is nearly temperature independent, so uncertainties in the temperature of 
ionized regions in galaxies are negligible. The Ha emissivity is parameterized directly in units of 
the luminosity density suggested by Gallego et al. (1995), but the uncertainty discussed above is 
included in the numerical coefficient. 

The uncertainties in the stellar emissivity (assumption 2) are summarized in Fig. 5. The 
error bars in the points from the Connolly et al. data are formulated in a different way from our 
evaluation of the AGN uncertainty. These authors corrected for survey incompleteness using a 
Schechter function with three different power laws to extrapolate to low luminosity. We show a 
rough la range in the emissivity based upon this estimate of the uncertainties. Our emissivity 
depends upon a fit to the data points in Fig. 5, and has small variations with the assumed 
cosmology. This introduces a small uncertainty of 0.01 dex in the calculation of Iq. 

The uncertainty inherent in our choice of opacity model, assumption (3) above, appears 
primarily in the degree of IGM attenuation of ionizing photons. If we focus only on the stellar 
contribution to the ionizing background, the number of He II ionizing photons produced is 
negligible. As a result, the contribution of diffuse ionizing radiation from the absorbers may be 
neglected. (Because H I recombination radiation is closely confined to the ionization edge, this 
diffuse radiation is quickly redshifted below vq.) Using Models 1 and 2 for the opacity instead of 
our standard model, we calculate slightly higher levels of Iq, by factors of 0.0068 and 0.0329 dex. 
A sample standard deviation of the models is ±0.017 dex and is probably an adequate assessment 
of the la uncertainty in Iq due to the opacity. If there is no opacity at all, To increases by 0.26 
dex. If the column density distribution for the absorbers in our standard model has an unbroken 
power law, with (3 = 1.3 rather than 1.5, the reduction in Iq is 0.21 dex. 

Our assumption (4), that the stellar ionizing radiation has a spectrum jy- 1 - 9 ^ - 2 (Sutherland 
& Shull 1999), yields a relatively small source of uncertainty. As a s is varied between 1.7 — 2.1, the 
specific intensity increases by factors of 0.045 to —0.043 dex. We therefore assign an uncertainty 
of ±0.045 dex. 

As shown in Fig. 8b, the contribution from stars at z > 2 to the present ionizing radiation 
background just above the Lyman edge is less than 3%. This assumes that the emissivity peaked 
at z ~ 2 and falls off at high redshift. If, as suggested by Pettini et al. (1998) and Steidel et al. 
(1999), there is little or no falloff in the star formation rate out to z ~ 4, the mean intensity at 
z = would increase by only 10%. This is a small effect because of redshifting and the fact that 
stars emit few photons more energetic than 4 rydbergs. Thus, stars at z > 3 would make little 
contribution to Iq at z = 0. Also, because stellar radiation does not doubly ionize He, diffuse He II 
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Lya and two-photon emission from absorbers make little contribution to Iq at z = 0. 

4. CONCLUSIONS 

In this paper, we have endeavored to make accurate estimates of the low-redshift intensity of 
ionizing radiation, arising from QSOs, Seyfert galaxies, and starburst galaxies. In performing this 
calculation, we found that we require accurate values of AGN emissivity and IGM opacity out to 
substantial redshifts. In other words, this is a global problem. 

Our new estimates of the ionizing emissivities of Seyfert galaxies and low-redshift QSOs 
were constructed by extrapolating ultraviolet fluxes from IUE to the Lyman limit. For starburst 
galaxies, we used recent Ha surveys together with an educated guess for the escaping fraction 
of ionizing radiation. The IGM opacity was derived from HST surveys of the low-redshift Lya 
absorbers and a new opacity model from Keck high-resolution spectra of high-redshift QSOs. By 
incorporating these ingredients into a cosmological radiative transfer code, we find that the ionizing 
intensity at z ~ has approximately equal contributions from AGN and starburst galaxies: 

/AGN = ^ 3 x 1Q -23 erg cm -2 g -i Hz -i gr -i) 

lf &r = (1.1 x 1(T 23 erg cm" 2 s" 1 Hz" 1 sr" 1 ) 

Taking into account uncertainties in the various parameters of the model (L mm , cx s , QSO 
luminosity function, Ha-determined star-formation history) we estimate uncertainties in these 
coefficients of 1.3^Q5 (for AGN) and l.lig 7 (f° r starbursts). For a spectral index a s ~ 1.8 (from 1 
- 4 Ryd), these values of Iq each correspond to one-sided ionizing fluxes <I>ion ~ 3000 photons cm~ 2 
s _1 . Allowing for statistical uncertainties, the total ionizing photon flux at low redshift probably 
lies in the range $i on = 2000 — 10,000 photons cm~ 2 s _1 , which is consistent with a number of 
recent estimates and measurements (see Table 1). 

The redshift evolution of the hydrogen photoionization rate, Tn\{z) is shown in Figure 9 for 
three cases: AGN only, starbursts only, and combined (AGN plus starbursts). These rates were 
computed using our standard emissivity models for AGN and starburst galaxies, except that the 
starburst emissivities were held constant at z > 1.7 to simulate recent measurements at high z. 
The starburst emissivities also assume (/ eS c) = 0.05. Although it is beyond the scope of this paper, 
the potentially dominant role of starburst galaxies in photoionization at z > 4 is apparent. 

Because the values in eqs. (14) and (15) include unacceptably large range for such an 
important quantity, it is worth discussing what might be done to improve the situation, both 
thoretically and observationally. Advances need to be made in the characterization of both 
ionizing emissivities and IGM opacities. The greatest uncertainty in the opacity model occurs in 
the range log Nhi = 15 — 17, where line saturation and small-number statistics make the Lya 
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surveys inaccurate. The imminent launch of the FUSE satellite will open up the far-UV band 
(920-1180 A) that contains Ly/3 and higher Lyman series lines. A survey of Ly/3 lines should 
allow better determinations of line saturation (doppler b- values). The FUSE spectra of AGN will 
also provide more accurate values of the flux near the Lyman limit and reduce the uncertainties 
introduced by extrapolating from spectral regions longward of 1200 A, as we have done with IUE 
data. 

To address the general problem of small-number statistics in the Lya absorbers, the HST 
Cosmic Origins Spectrograph (Morse et al. 1998), scheduled for installation on HST in 2003, 
should be used for a QSO absorption-line key project. A full discussion of the advantages of this 
project is given in Appendix 1 of the UV-Optical Working Group White Paper (Shull et al. 1999). 
Because COS will have about 20 times the far-UV throughput of the previous HST spectrographs, 
GHRS and STIS, and offers velocity resolution 10 times better than that of FOS, it will provide 
much better statistics on the distribution of H I absorbers in both space and column density. 
The most useful COS surveys will be of low-redshift Lya lines, particularly the rare "partial 
Lyman-limit systems" (16.0 < log Nhi < 17.5). Accurately characterizing the distribution in 
column density and the evolution in redshift of these absorbers would remove a large part of the 
uncertainty in the IGM opacity model. 

The emissivity models for both AGN and starburst galaxies also need improvement. Although 
we have used current surveys of Seyfert galaxies and QSOs, we may have missed certain classes 
of sources that are strong emitters in the EUV. We believe that BL Lac objects contribute less 
than 10% of Seyferts, based on estimates of their luminosity function and space density. On the 
other hand, Edelson et al. (1999) suggest that narrow- line Seyfert 1 galaxies may account for 
~ 50% of the EUV volume emissivity in the ROSAT Wide-Field Camera sample. It is not clear 
whether these Seyferts are captured in the Cheng et al. (1985) luminosity function, but their 
ionizing spectra might be higher than that derived from an extrapolation of their UV fluxes. For 
low-redshift starbursts, two recent surveys (Gallego et al. 1995; Gronwall 1998) derive comparable 
values of Ha luminosity density, although even the latter (KISS) Ha survey may still be incomplete 
at the faint end. At higher redshifts, the QSO luminosity function is uncertain, owing to the effects 
of dust (Fall & Pei 1993; Pei 1995) and faint-end survey incompleteness. QSO surveys by GALEX 
(Galaxy Explorer) in the ultraviolet and by the Sloan Digital Sky Survey in the optical may clarify 
the AGN luminosity functions. However, it is worth repeating that, owing to redshifting and IGM 
opacity, AGN and starbursts at z > 3 contribute less than 10% to the value of Iq at z = 0. 

For an accurate emissivity density, what is needed most are surveys at z < 1 of AGN and 
starburst galaxies. Even after we ascertain the space density and Ha luminosity functions of 
star-forming galaxies, we still need an accurate measurement of (/ csc ), the fraction of stellar 
ionizing photons that escape the galactic H I layers into the halo and IGM. Here, we have relied 
on recent theoretical work (Dove et al. 1999) and Ha observations of gas in the Magellanic 
Stream (Bland-Hawthorn & Maloney 1999) that suggest (/esc) ~ 0.03 — 0.06. However, access 
with FUSE to the far-UV spectrum at 920-950 A allows a direct measurement of the escaping 
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EUV continuum from starbursts at redshifts z « 0.05. This work will extend the studies with the 
Hopkins Ultraviolet Telescope (HUT) of leaky starbursts (Leitherer et al. 1995; Hurwitz et al. 
1997). 

Finally, we eagerly await new measurements of the ionizing photon flux, <E>i on via several 
direct and indirect techniques. These methods include improved Fabry- Perot measurements of 
Ha from Galactic high- velocity clouds (Tufte et al. 1998; Bland-Hawthorn & Maloney 1999), 
and UV absorption-line measurements of ionization ratios such as Fe I/Fe II and Mg I/Mg II 
that constrain the far-UV radiation in the 0.6-1.0 Ryd band (Stocke et al. 1991; Tumlinson et 
al. 1999). A new absorption-line key project with HST/COS could make precise estimates of Iq 
from the proximity effect. As a result of the new surveys and ventures mentioned above, it should 
be possible, within five years, to determine the local ionizing background to < 30%. With good 
fortune, these measurements and the theoretical models will agree to a level better than described 
here in Table 1. 
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Table 1 

Measurements and Limits of low- z Ionizing Background 1 



Technique 



$ ion (cm 



2 s- 1 ) 



Reference 



Ha Fabry-Perot 
Ha Filter Images 
Ha Filter Images 
Ha Filter Images 
H I Disk Edges 
H I Disk Edges 
Prox. Eff. (z) = 0.5 



< 3 x 10 4 

< 1.1 x 10 4 

< 8.4 x 10 4 

< 9 x 10 4 
(0.5 - 5) x 10 4 

(1 - 5) x 10 4 
(0.05 - 1.0) x 10 4 



Vogel et al. (1995) 
Donahue et al. (1995) 
Stocke et al. (1991) 
Kutyrev & Reynolds (1989) 
Maloney (1993) 
Dove & Shull (1994a) 
Kulkarni & Fall (1993) 



1 f&ion is the one-sided, normally incident photon flux in the metagalactic background, related to 
the specific intensity at Lyman limit by $i on = (2630 cm~ 2 s _1 )/_23(l-8/a s ) - see eq. (1) in text. 



Table 2 
Low-z Opacity Models 



Model 


iV min (cm 2 ) 


AW (cm 2 ) 


A l 


ft 


7* 


Ci 


Si 


Standard 




10 14 


0.105 


1.63 


0.15 


0.010 


-2.81 




10 14 


10 17 


0.501 


1.50 


0.15 


0.553 


-2.73 




10 17 


1Q 22 


0.159 


1.50 


1.5 


0.263 


-1.04 


Model 1 




10 15.5 


0.105 


1.63 


0.15 


0.048 


-2.81 




10 15.5 


1Q 22 


0.159 


1.50 


1.5 


0.411 


-1.38 


Model 2 




10 15.5 


0.105 


1.63 


0.15 


0.048 


-2.81 




10 17 


1Q 22 


0.159 


1.50 


1.5 


0.263 


-1.04 
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Fig. 1. — Evolution of the ionizing background with redshift, assuming the modified Pei luminosity 
function for AGN described in § 2.1 (0,q = 0.2, h = 0.5, auv = 0.86, a s = 1.8). The opacity is a 
hybrid model consisting of Model A2 from FGS for z > 1.9, and our standard model for the low 
redshift opacity for z < 1.9. This calculation incorporates the full radiative transfer and reradiation 
of ionizing photons by absorbers described in FGS. Top panel: 3 < z < 5. Bottom panel: < z < 2. 
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Fig. 2. — Correlation between L912 and Lb for the 27 Seyferts in our subsample. A least-squares fit 
gives a slope of 1.114 ± 0.081 with \ 2 = 65.5. For comparison, the dashed curve shows a constant 
optical spectral index of auv = 0.86 between 4400 A and 912 A. 
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Fig. 3. — AGN luminosity functions. Circles show the Seyfert galaxy sample of Cheng et al. (1995) 
from which the 27 Seyfert galaxies in common with our IUE-AGN database (Penton et al. 1998) 
were taken (see Fig. 2). Squares show the Seyfert galaxy sample of Kohler et al. (1997), and dotted 
curve represents a power law fit $(£#) oc L~^ 2 ' 2 to these latter results. The solid and dashed curves 
are an extrapolation to z = of the Pei (1995) luminosity function for h = 0.5, 0,q = 0.2 and 
h = 0.5, Qq = 1, respectively. Note that at z = 0, the AGN fitted by Pei only went down to 0.7L* 
corresponding to Mb = —23.5. 
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Fig. 4. — Source emissivity due to AGN at redshift z for modified Pei open model (f2o = 0.2, 
h = 0.5, and auv = 0.86) (solid curve) and modified Pei closed model (Qo = 1, h = 0.5, and 
auv = 0.86; correction for open universe line element) (dashed curve). Dotted curves frame our 
la range in emissivity for the open model, including corrections for completeness of the luminosity 
function, and a variation of the L^/Lb ratio by one standard deviation in spectral slope and 
dependence on Lb- 
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Fig. 5. — (a) Source emissivity due to galaxies at redshift z, tied to Gallego et al. (1995) and 
Connolly et al. (1997). Assumes ft = 0.2, h = 0.5, and {f esc ) = 0.05, T = 2 x 10 4 K, {hv) = 22 
eV. Dotted curves frame our la range in emissivity. (b) Same as (a), except Qq = 1 
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Fig. 6. — Differential continuum opacity, dr e g/dz, at v = uq and z = 0, versus the assumed upper 
limit on the column density distribution. We consider HST/FOS sample 5 of Weymann et al. (1998) 
of 465 Lya absorbers with W\ > 240 mA, for which dj\f jdz = 30.7 ± 4.2. We model these lines as 
a distribution, N^j, of column densities from Nhi = 10 14 cm~ 2 up to Nhi, u - I n ascending order, 
these curves assume (3 = 1.7, 1.5, and 1.3. 
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Fig. 7. — Differential continuum opacity, dr e s/dz, at v = uq versus redshift z for models for the 
column density distribution of Lya forest and Lyman limit systems. Solid curve is our "standard 
model" described in Table 2. Dotted curve is our Model 1, and short-dashed curve is our model 2. 
Dot-dashed curve is the extrapolation to low z of Model A2 from FGS. Long-dashed curve is the 
opacity assumed by Haardt & Madau (1996). 
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Fig. 8. — (a) Cumulative value of Iq{z = 0) arising from AGN at redshifts greater than z. These 
intensties are calculated using the approximate method (eqs. (3-5]) that neglects diffuse ionizing 
radiation from the absorbers. Solid curve assumes f2o = 0.2, h = 0.5; dotted curve assumes 
f2o = 1, h = 0.5. Standard model for opacity is used for both. At z = 0, Iq = 1.09 x 10~ 23 
and Iq = 1.13 x 10 -23 for these respective cases, (b) Cumulative value of Iq arising from galaxies 
at redshift z, assuming fi = 0.2, h = 0.5, (solid curve) or fl = l,h = 0.5, (dotted curve), 
(fesc) = 0.05, and standard model for opacity. At z = 0, I = 1.13 x 10~ 23 and I = 1.07 x 10~ 23 
for the respective cases. 
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Fig. 9. — Hydrogen photoionization rate versus redshift for three cases: AGN-only (dashed curve), 
starburst galaxies only (dotted), and combined (solid curve). Each case assumes the standard 
models for emissivity (§ 2.1, 2.2) and opacity (§2.3 and Fig. 1 caption). The AGN emissivities 
are the same as solid curve in Fig. 4, while starburst emissivities are same as solid curve in Fig. 
5a, except that the comoving emissivity is held constant for 1.7 < z < 4.0 to simulate recent 
measurements of star-formation history. 



